Dynamical Vortices in Superfluid Films 
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The coupling of superfluid film to a moving vortex is a gauge coupling entirely dictated by topology. 
From the definition of linking number, one can define a gauge field _4 M , whose (2+l)-dimensional 
curl is the vortex 3-current J M , and to which the superfluid is minimally coupled. We compute 
the superfluid density and current response to a moving vortex. Exploiting the analogy to (2 + 
f )-dimensional electrodynamics, we compute the effective vortex mass M(u>) and find that it is 
logarithmically divergent in the u — > limit, with a constant imaginary part, yielding a super- 
Ohmic dissipation in the presence of an oscillating superflow. Numerical integration of the nonlinear 
Schrodinger equation supports these conclusions. The interaction of vortices with impurities coupling 
to the density also is discussed. 

PACS numbers: 67.40.V 
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I. Introduction 

In this paper we investigate the effective action and dy- 
namics of vortices in compressible superfluid films at zero 
temperature. In an incompressible two-dimensional (2D) 
superfluid, vortices behave as massless charges in a uni- 
form magnetic field - their motion is along an equipo- 
tential, the sum of logarithmic contributions from each 
of the point "charges" (I) . The Lagrangian for a charge- 
neutral system of vortices may be written 



L = — np^^^ rii XiYi + ^ — riiUj In \Xi — Xj 



2tt 



where k = h/m is the rotational quantum, p is the 
bulk density (or superfluid density Q), rij is the inte- 
ger "charge" and Xi the position of the i th vortex. The 
equations of motion, 



— * fC x — * 

Xi = — > n, 

j 



zx (Xj-Xj) 

\Xi-Xj\ 2 



preserve the total potential energy of the vortices, which 
of course is just the kinetic energy of the superfluid. 
These equations are first order in time - there is no iner- 

tial term £\ \MiX\ in L. 

In a compressible superfluid, the speed of sound c is finite. 
This leads to retardation effects in the vortex dynamics. 
Furthermore, accelerating vortices may radiate phonons, 
leading to dissipation. Both effects are described by a 
complex frequency-dependent mass term M(uj), derived 
below. This physics is present in granular films and 
Josephson junction arrays as well J3|. The basic idea is 
to integrate out the phonons, which represent a bosonic 
bath to which vortices are coupled, in the spirit of Ref. 
Q , and thereby derive an effective action for the vortices 
alone. 

This paper is organized as follows: In section II we derive 
the analog of backflow for moving vortices in superfluid 



films. In section III we review the correspondence be- 
tween superfluid dynamics and electrodynamics in two 
space dimensions and show how the results of section II 
may be obtained by a Lorentz transformation of a static 
vortex solution. Self-interaction effects, vortex mass, and 
dissipation are discussed in section IV. In section V, we 
report on the results of numerical simulations of a vortex 
in an oscillating superflow, from which we can extract 
M(oj) and compare with theoretical predictions. Section 
VI discusses the interaction of vortices and dynamical 
impurities. 

II. Analog of Backflow for Dynamical Vortices 

Consider a vortex moving in a (2 + l)-dimensional Bose 
fluid. The only information we have about the vortex is 
that it is a point object which accrues a geometric phase 
of 27r in the many body boson propagator each time it 
encircles a boson. We write the vortex current density as 



CK 



Mr)) 



where r parameterizes the vortex "world lines" X, (r), 
which are one-dimensional filaments running through 
(2+l)-dimensional spacetime. The many boson La- 
grangian is written 



£ = 



•Kj I ) ~f" -^top 



i<j 



where we assume a simple generic interacting Bose fluid 
(isotropic, single component). Here, Ct op is the topolog- 
ical term in the Lagrangian which counts the winding 
number of the vortices relative to the bosons. This is 
explicitly written in terms of the linking number of their 
trajectories 0, 



Stop = Jd 2 xdt C t0 p = 2nh N iink 



d 2 xdtj>*^^J X = jil-rdi : (1) 
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where the boson mass current density is j M = (cp,j), 
and where d v jd 2 is a formal expression for a nonlocal 
operator. Vortex current conservation, d^J^ — 0, allows 
one to construct a gauge potential A^ whose curl is J M 



J^ = -ce„ vX d v A x , (2) 

and thereby express the linking number as a local in- 
teraction between the boson current density and the 
vortex gauge potential A^ . 

The time-dependent Hamiltonian for the bosons in the 
presence of moving vortices is thus 



H{A») = H(0) - / d 2 x3l{x)A»{x,t) 



+i d 2 xp(x)A 2 (x,t) (3) 



where p is the boson density and j p is given by 



i p 
JO 



(x) = cp(x) = mcS^S(x — xi) 

i 



The gauge invariant boson current density is then writ- 
ten H: 
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At this point we have reproduced the well-known analogy 
between a rotating superfluid and a superconductor in a 
magnetic field ||] . What is new here is the explicit gauge- 
covariant, time- dependent description, through Eqs. |^ 
and [|, of the coupling of a superfluid to vortices, which 
are quanta of rotation. 

The linear response of the boson system is given by the 
Kubo formula, 



where 



d 2 x' dt' K^(x, t; x", t') A v (tf, t') 



K^(x,t;x',t') = jr(\j»(x',t),j»(x / ,t')}) G(t-t') 

-{p(x, t)) <f M „ (1 - Sua) S(x - x') S(t - t') . 

The spatial part of K^ v may be written in terms of longi- 
tudinal and transverse components in Fourier space, viz. 

K ij (q) = ~q l ?K^{q) - {g* + q^)K \{q) . 

Gauge invariance, K liu (q)q u = 0, may be used to relate 
the 00 and Oi components to K ii : 



K 00 (q) 



c 2 \q\ 2 



K l{ (q) , K i0 (q) = -^-K l{ ( q ). 



At zero temperature, the single mode approximation 
(SMA) gives for the density response JToJ 



K s r(q) = ! 4 l s(q: 
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to + c\q\ + ie lo — c\q\ + ie 



where p — (p) is the average mass density and S(q) is 
the ground state static structure function. Note that 
K^ MA (q) — p, since the phonon is purely longitudinal. 
Recall that lim^o h\q\/2mcS(q) = 1 (To). 
Static Vortex - We choose a gauge in which A = 
and V x A = nK,5(x)z, which is satisfied by A(x,t) — 
nnzxx/2Tr\x\ 2 . Now A{q,uS) = — inn (zxq/\q \ 2 )-2ttS(uj) 
is purely transverse, so the density response vanishes and 
the current density response gives the usual (j {1) (x)) — 
nnpz, x x/2ir\x\ 2 . The absence of a density variation 
in response to the vortex seems to contradict the result 
that 5p(r)/p = —n 2 Ksph 2 /2m 2 \x\ 2 far from a vortex 
of strength n, where K$ is the adiabatic compressibility. 
However, the n 2 dependence tells us that this is a non- 
linear response. The second order response is formally 
written 



d 2 x'dt' \ d 2 x"dt" R al 3 7 (x; x'; x'^A^x^A^ix"). 



may be divided into two contri- 
butions. The first, R 1 ^, is the second order nonlinear 
susceptibility arising from the linear j^A^ coupling in 
H. The second, R 1 ^, is the linear susceptibility arising 

from the \pA 2 term - this is given by the density-density 
correlation function, so 

R^p^x; x'\ x") = --^K o(x, t;x',t')5(x' - x") 
x5(t'-t") 6^(1-600) . 

It is easy to check that the nonlinear response arising 
from i2jjg 7 exactly reproduces the asymptotic density 
variation due to the vortex, and that does not con- 

tribute to 0(l/\x\ 2 ). 

Moving Vortex - Consider now a vortex moving with uni- 
form velocity: X = vt. We choose the gauge ^4° = and 



A(q, uj) = —inn 



z x q 



■ z • v x q 



■2ir 6(lu — q-v) 



Note that this is no longer purely transverse, so there will 
be a linear response of the density to the moving vortex: 



= 8inn 2 pS(qh 



z ■ v x q 



8(uj — q-v) . 



\q\ {cqf-{v.qf 
Further assuming \v/c\ <C 1, we obtain, at large distances 

nnp z • v x x 



(p (1 W)> 



2ttc 2 \x-vt\ 2 ' 



(4) 
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which is identical to the result obtained by Duan [12]. 
The current density in the presence of a moving vortex 
is similarly computed and found to be 



, . . . nnp z x R 



1 - 



2c 2 



(V ■ Rf - (v x Rf 



valid to order ?7 2 /c 2 , with R = x — vt. 
Contrast with Impurity BackRow 

It is important to contrast this behavior with the stan- 
dard picture of backflow in neutral systems. A local per- 
turbation coupling to the superfluid as 

H' = [ d 2 x p(x) U(x, t) 



leads to a superfluid response 



1 



(j^(q,co))=—K M (q,oj)U(q,Lo)e l . 
c 

Using the SMA response functions, and assuming 
U(x,t) = Uq5(x — vt), one obtains 



( P ^(x,t)) = -^S(R) 
pU v 



0^x,t)) = -^f 2 {v-2( v .R)R} 



(5) 
(G) 



The density response is purely local, in contrast to that of 
Eq. ^, and the current, which vanishes in the static case 
v = 0, is dipolar and falls off as l/R 2 . The superfluid- 
vortex gauge coupling leads to a much different linear 
response. 

III. Analogy to QED 2+ i 

These results may be understood in terms of the well- 
known correspondence between (2 + l)-dimensional su- 
perfluids and quantum electrodynamics fl3|-|l7|l, which 
we now review. One starts with the standard Ginzburg- 
Landau Lagrangian density in the presence of an external 
gauge field 



C[W*, ^1 = &*{ihd t + eZ°)<P- — 

2m 



i c 



(7) 



For a superconductor, Z^ would represent the electro- 
magnetic gauge potential, while in our case it can be 
used to describe an externally imposed current. At this 
point, the "charge" e and velocity c are arbitrary param- 
eters; we will take c to be th e spee d of sound, defined 
below. One substitutes \P = \fp~fm e ie e 1 ^, where 9{x,t) 
is a smooth "spin- wave" field and X(x,t) the singular 
vortex field, which satisfies J^(x) — (hc/m) e^ uX d v d x X. 
This gives 



£' = -HE (d t e + d t x - |z ) - ^4 f V0 + vx + 2 

m V n / 2m z V he 

-2 



H (Vp) 2 - ±(p-p? 



8m 2 p 



m- 



after subtracting a time derivative term. Decoupling the 
V# + VX + eZ/hc) term, one arrives at 

c „ = _Wf^ 9 + ^ x+ eg)_*P ( d9 + dx _!l Z o 
m V Tic / m V Ti 

Q 2 h 2 (Vp) 2 A 



2p 8m 2 p 



Integrating over the spin wave field 6{x, t) now generates 
the constraint V • Q + d t p = - Q is the mass current 
- which is satisfied by introducing the gauge field A 11 = 
(A , A), where 

(p - P, Q) = -- (z • V x A, cz x V4° + z x d t A) . 
c 

The coupling of the vortex current to this gauge field is 
due to the term 



h 

m 



(Q ■ VX + pd t x) = ^ e^ A (A„ + afl )d v d x X + d( ■ ) 



m 
1 



= -pJ» (A^ + a^ + di ■ ) 



where the nondynamical gauge field generates a static 
magnetic field, i.e. 

e = -Va°---^=0 , 6 = z-Vxa=-c, 
c at 

which is satisfied by the gauge choice a = \c(y, — x), 
a = 0, for example. The coupling between the gauge 
fields A^ and Z v is then given by 

- (pZ Q - 1 -Q.Z]=^- e^Z^Ax + o A ) . 
m \ c ) mc 

Finally, one introduces the field strength tensor 

/ E x Ey\ 

Fp, = dpAv - d v A^ =[-E x —B , 
\-E y B / 

and obtains 

C cS = -| V" - p - J»(A„ + a„) 

+^-e? vX Z ll d v (A x +a x ) 
mc 



+ 



p BE 2 p (£VB) 2 



2c(l-B/c) 8 (1-B/c) 



(8) 



where c = y // 2\p/m 2 is the speed of sound and ^ = Ti/mc 
is the coherence length. Note K° = m 2 /2p 2 X is the bare 
compressibility. 



3 



The Lagrangian density C c s describes "charged" parti- 
cles (vortices) moving in a background "magnetic field" 
— cz, minimally coupled to a dynamical gauge field 
(note that this is not the gauge field defined in Eq. ||). 
That the background magnetic field is the average bo- 
son density is of course due to the fact that the vortices 
see the bosons as sources of geometric phase. This was 
recognized by Haldane and Wu |l8|], who computed the 
Berry phase accrued by a vortex as it executes adiabatic 
transport in the superfluid film. If the vortex position is 
£, and the adiabatic wavefunction is | ^ ) , then 



7c 



>|v f i*; 



-2tt^- S, 



C ' 



where S c is the area enclosed by the path C along which 
the vortex travels. Note that this immediately tells us 
that vortices experience a Lorentz force when moving 
through a superfluid film. In the nonrelativistic limit, 
their dynamics can be described in terms of vortices being 
advected in each other's flow field, or as charged particles 
in a background magnetic field moving under the influ- 
ence of each other's electric field. A vortex-antivortex 
pair, for example, behaves like an exciton in a magnetic 
field. 

The linearized, long wavelength Lagrangian density is ob- 
tained by dropping terms in C e s which are higher than 
second order in the field strength or which involve higher 
derivatives acting on the field. One is then left with 
(2 + l)-dimensional quantum electrodynamics, 



-QED 



4 " 



c Vm / 



in the presence of a uniform background magnetic field. 
Note that the gauge field A M couples to a sum of the 
(quantized) vortex current density J' 1 and the (not quan- 
tized) externally imposed "current" (e/m)e^ x d tl Zx, 
which could represent a global rotation of the system |j| . 
When no vortices are present ( J M = 0), one can integrate 
out the gauge field A^ to obtain 



p / e \ 2 fcPkcLu 
= 2 \mc) J~(2nf 



k^k 1 



Z„(k) 



m 



where k 2 = k^k u 



k 2 , and setting e/rac = 1 



one can read off the response tensor 



K^{k) = p 



k»k v 
k 2 



When Z M = 0, one has a theory of vortices minimally 
coupled to the gauge field A^, and the action extrem- 
izing equations for the fields are Maxwell's equations: 
dpF*" = J v /c, or 



V-E = -J° 

c 

- - IdB 

z-VxE = — 

c at 

VB x z = -J-\ — . 

c c at 

(Note V- B = dB/dz = 0, trivially.) For the vortices, 
one has the Lorentz force law, 



dX l 
dXf 



x E[Xj) 
-B(X t ) 



(9) 



which says that the vortices move perpendicular to the 
local electric field, with a magnetic field strength of c — B 
which is the sum of a uniform background contribution 
(the average boson density) and a dynamical contribution 
(due to fluctuations in the boson density). 

Lorentz Transformations 

To investigate the effects of moving vortices, it is useful to 
appeal to the Lorentz invariance of £qed and transform 
static solutions pl| . Recall that in (2 + l)-dimensions, 
the Lorentz group has three generators, corresponding to 
two boosts and one rotation. The general boost trans- 
formation is written 



( 7 

\7A/ 



^B 2 



^B 2 



1 



(10) 



with 8 = v/c and 7 = 1/y 1 — B 2 . Applying the Lorentz 
transformation z'^ — L p ' v z v to the coordinates and field 
strength tensor gives the familiar results 



x = "fX + 7/3 • x 
7-1 



7x u /3 



P 2 



(B ■ x)B + x 



and 



E' = jE — 
B' 



7' 



1 



~-{B ■ E)B + ~fBi x 3 



B 2 

7z • (3 x E + -fB . 
We may now transform solutions 

{x^,J",F^} — ► {x'>*, J^,F'^} . 

A static charge 1 vortex generates an electric field E = 
kx/2tt\x\ 2 and a magnetic field B = 0. Upon applying 
the boost of Eq. [l^, we obtain (dropping primes), 



E = 



B 



ry K x ± z x j3 + (x,, - Bx°)j3 

2^ x\ +^f{x {l -Bx ) 2 

7« 0x ± 
^xl+^ix^ - Bx*) 2 



(11) 
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where we have written X* = x,,@ + x ± z x 0. 

Now the rules for translating from E and B to j and p 
are as follows: 

J=pi,xE p = p(l- B/c) 

We now see that the linear response formulae (Eq. || and 
accompanying discussion) exactly reproduce these results 
to lowest order in (3. The vortex velocity is the the ratio 
j/p, which is the content of the Lorentz force law, Eq. |j. 

A Tale of Two Vortices 

Consider now an elementary vortex-antivortex pair sep- 
arated by a distance a. We choose X±(t) = vtei ± iaeV 
Computing the electric and magnetic fields at one of the 
singularities due to the presence of the other is easily 
accomplished with the Lorentz transformation. One ob- 
tains 



E(X i 



"fK 

2na 



B(X d 



J3jk 
2na 



But if in the moving frame the vortices are stationary, we 
must have that c(3 = dX±/dt, which leads to the result 



/3(a) 



^ 2 



1/2 



Thus, at large separations the pair's velocity is c£/a, 
but at smaller separations the velocity asymptotically ap- 
proaches the sound speed c. We stress that this is true for 
the model defined by Cqed, where the vortices have no 
core. The naive expression v(a) = c£/a begins to break 
down at distances on the order of £, where a proper ac- 
counting of the terms neglected in the QED action must 
be taken in order to reproduce the correct core structure. 

SuperQow and the Magnus Force 

In the above examples, we derived results for a moving 
vortex in a stationary superfluid. In this section, we make 
a Galilean transformation (the original theory is Galilean 
invariant!) in order to discuss what happens to vortices 
in the presence of a background superflow. Starting with 
the Galilean-transformed Lagrangian density 



C = ih&*d f &+ihv-&*V&-— |WT-A ( l*f - — 

m 



2m 



and proceeding as before, one derives the effective La- 
grangian density 



(E - Bz x (3) 2 



1 



B/c 
1 



B 2 



P (gva) a 

8 1 - B/c 



pJ^A, 



where now generates a static electric field as well as a 
static magnetic field: 



e — — cz x (3 , b = —c . 

The Lorentz force due to e is the Magnus force. The 
vortex equation of motion is 



dX l 
dXf 



When E{X l ) + 
it is stationary. 



= + 



i x E{Xj) - @B{X l 
c-B(X) 



0, the forces on vortex I cancel, and 



IV. Self-Interaction, Inertial Mass, and Dissipation 

The analogy to electrodynamics suggests that there 
should be an electrodynamic contribution to the mass 
and retardation effects, as there are in (3+l)-dimensional 
classical electrodynamics p2[ . In the superfluid, this is 
due to the phonon cloud carried by the vortex - a po- 
laronic effect. However, as we've seen, the coupling of 
vorticity to superfluid density and current fluctuations is 
a gauge coupling which is rather different from the local 
density coupling used in conventional polaron theories. 
Still, this coupling is of the general form considered in 
Ref. [|J, i.e. an external coordinate (the vortex position) 
coupled to a bath of oscillators (the phonons) . 
We wish to integrate out the dynamical field A^ corre- 
sponding to the phonon degrees of freedom and obtain 
an effective action for the vortices. Working in Lorentz 
gauge (d^A^ = 0), we integrate out the Gaussian field 
A^ in £qed by solving the equations of motion, yielding 

A^ = - W 1 ^ 

c 

and an effective action for the vortices of 



5, 



eff 



d 3 x dV J" (a;) □~ 1 (ar, x') J^x') 



P 



/ d 6 x J»{x)a^{x) . 

The inverse D'Alambertian has the form 

1 0{x° -x'° ~\x-x'\) 



□ _1 (x,x') 



27T ^( x _ 3/0)2 _ |£_£'|2 

in (2 + l)-dimensions and 

5{x° -x /0 -\x~x'\) 



LT" 1 (a? ) af') = 



4-7T | X — X 1 1 



in (3 + l)-dimensions. Thus, in contrast to the case of 
three spatial dimensions, where □ vanishes unless x—x' 
is light-like, in our (2+l)-dimensional cased -1 is nonzero 
everywhere inside the light cone [^9|. The finite sound 
speed c leads to retardation effects. One might naively 
think that this would lead to the collapse of the vortex- 
antivortex pair, since the vortex should 'see' the antivor- 
tex at earlier times and vice versa. However, although the 
potentials are retarded, the fields of a uniformly moving 



5 



charge point to the instantaneous position of the charge 
(as we've derived above), and so for the special case of 
a uniformly moving vortex-antivortex pair, there is no 
apparent time delay pofl . 

The self-interaction part of S c g for a vortex of strength 
n is 



Sscli — — 



2 2 — 

n k p 

47TC 



du I da 

oo Jo 



1 - (3{u) ■ (3{u + a) 
6(o-- \X(u + a)-X(u)\) 



\X(u + a)-X(u)\* 



(12) 



where we've taken x° (u) = uasa parameterization of the 
vortex world line, and (3{u) = dX{u) / du. If the integrand 
in S'seif were well-behaved and allowed an expansion of 
the form 



= / dt {-m c 2 + imiu 2 (£) + . . .} 



>solf 



then we would associate a rest mass with mo and an 
inertial mass with mi. The remaining terms, involving 
higher derivatives and powers of the velocity, would be 
negligible in the Newtonian limit. However, the integral 
over a in Eq. [l2] diverges logarithmically, both for large 
and small a. The small a divergence is remedied by a 
proper treatment of the core structure, which lies be- 
yond the QED approximation. The large a divergence, 
on the other hand, is real. In the case of the parameter 
mo, this is to be expected, since we know the energy of a 
static vortex diverges logarithmically with the size of the 
system, owing to the slow fall-off of the current density 
\j\ oc 1/r. In the electrodynamic language, the energy 
density, £ = \p (E 2 + B 2 ), dies off as 1/r 2 in the vicinity 
of a static vortex, yielding a logarithmic divergence when 
integrated over the system. In (3 + l)-dimensional classi- 
cal electrodynamics, by contrast, the energy density dies 
off as 1/r 4 , and there is no infrared divergence. Now we 
ask whether mi is finite. The answer again is no. Since 
iS S eif is a Lorentz scalar, for constant (3 one has 



5 so if(/3) = -VI -fi 2 / dtmoc 



which says that mi = mo. As recently emphasized by 
Duan flljQ, this may be understood in terms of the 



density variation (p (1) 



(X z-vxx/ x- 



-vt\ 2 , which produces 



a logarithmically infinite energy shift 



AE 



1 



2K s p 2 



d 2 x [5 P {x)f 



Similarly, the total momentum P of the moving vortex, 



P (1) (t) =m \d 2 x 2 ^' l) 



p(x,t) 



Frequency-Dependent Inertial Mass 

In this section we compute the low frequency inertial 
mass of a single vortex and find that it is frequency- 
dependent and logarithmically divergent as u> — » 0. We 
start by expanding the self-interaction contribution 5 sc if 
for a single vortex [f24f: 



Ssoif = - . / du I da 



Aire 



oo JO 



1 - f3(u) ■ f3(u + a) 



1 , \X{u + a) - X(u) 

17 

oo 



2a 3 



= S,,.,u, +in 2 / > ^-M'(u J )oj 2 \X(u)\ 2 + 



where <S s t a tic is action for a static vortex of strength n. 
The quantity M'(uj) is found to be 



M'(w) = 2^ 



D ds 



1 — COSLUS 



COSCJS — 



(13) 



M ci(M<5) + 



1 — cos lo8 sinw(5 



j2 6 2 



where u = Tr£ 2 p = nTi 2 /2\ is the "core mass" of the 
vortex p2], and ci (z) is the cosine integral |E5[. We 
have introduced an ultraviolet temporal cutoff 5 « £/c 
to regularize the s integrals. This crudely accounts for 
the core structure of the vortex which lies beyond the 
approximation afforded by £qed- The important point 
is that the infrared divergence of (S— Static) is suppressed 
by the finite frequency oj, leading to a low frequency mass 
n 2 M'(bj) which diverges logarithmically as lo — * 0: 

M'M - p{-HW) -(C + |) + ^ 2 <5 2 + ...}, 

where C = 0.577215 ... is Euler's constant. 
The effective Lagrangian C e g does contain a term 
— hp (£Vi?) 2 which is dropped in the long wavelength ef- 
fective theory but is still quadratic in the field strengths. 
(It also breaks Lorentz invariance.) Retaining this term, 
there are no ultraviolet divergences, and the mass is 



M(u>) 



In 



yj A(tu) + 1 \ iirp sgnw 
M'{uj)+iM"{uj) 



(14) 



A(u) = 1 + 



uj 2 e 



(15) 



where we have now included the imaginary part M"(u)). 
The logarithmic divergence at small u> is still present, but 
in the large lo limit we find that M'(u>) vanishes as ui~ 2 
and M"{uj) as uo^ 1 . The Fourier transform of M(uo) is 
then causal: 



diverges logarithmically |12 



G 



M(t) 



[I (ct/0 - L (ct/0] 6(f) 



7T/iC 

^[i-^/o + i(c«/o a 



(t-0) 



[1 - (£/ci) 2 



(i -> oo) 



where Jo (2) and L (z) are modified Bessel and Struve 
functions, respectively p5fl . The logarithmic frequency 
dependence has previously been obtained by Eckern and 
Schmid ||, who investigated vortices in granular films, 
and by Stamp, Chudnovsky, and Barbara [^8j in the con- 
text of magnetic domain walls. 

When several vortices are present, the effective vortex 
action for low frequencies becomes (see the Appendix) 



' dui 
2/ 2ir 



J~]mX'i(u 



pn 

47T 



dt ^2 n i n 3 ln \Xi(t) - Xj(t)\ 



-kp J dt Yl 



mXiffiYiit) , (16) 



where the first term arises from an expansion of J' 1 □ 1 
in terms of the vortex coordinates themselves p9[] , and 
where the prime on the sum is a zero total vorticity re- 
striction: X)i n « = 0- Notice that the first ("kinetic") 
term, discussed in the Appendix, involves only the to- 
tal dipole moment operator D(t) — ^\ n>iXi(t); this fact 
is intimately connected with Galilean invariance and the 
stability of superflow in the absence of disorder. Con- 
sider, for example, an elementary (n = ±1) vortex- 
antivortex pair. Let X = \{X+ + XJ) be the "center 

of mass" (CM) coordinate and x = (X + — X-) be the 
dipole moment. The CM coordinate appears only in the 
Berry phase term of the Lagrangian, which is 

L B = -Kp{X + Y+ -X_y_) = -Kp{Xy-Yx) 

up to a total time derivative. Thus, a path integral over 
the CM coordinates generates a delta function at each 
time step, enforcing dx(t)/dt = always. 

Relative Importance of Inertial Terms 

To investigate the importance of inertial terms relative 
to those arising from the Lorentz force, we consider the 
response of an isolated n — 1 vortex to a time dependent 
field e(t), which might represent a sudden switching on 
of a superflow which will accelerate the vortex |1 or an 
oscillating superflow. We find that the velocity V{u>) = 
— iu>X(u>) satisfies 



V(u) 



1 



i x e(u>) 



where the dimensionless function 



ir(u>) 
1 -r 2 {io) 



e» (17) 



wM{u) _ uj£_ M(ui) 
Kp 2c p 



shown in Fig. 1, describes the inertial and frictional as- 
pects of the vortex's motion. 

We now see that inertial effects will be relatively unim- 
portant at frequencies where the "inertial parameter" 
r(u>) is small. Note that r'(u>) vanishes both for very 
low and very high frequencies, as shown in Fig. 1. Tak- 
ing c ~ 100 m/s and £ - 5A, one obtains a characteristic 
frequency too = c/£ ~ 10 11 Hz. At low frequencies, both 
real and imaginary components of r(u) are small, and 
inertial effects are relatively unimportant. From Eq. [Tt], 
we find that an elementary vortex in an oscillating su- 
perflow will move at a Hall angle 6h(uj) — tan -1 \r(ui)\ 
relative to the z x e direction. The power dissipation per 
unit frequency is given by 

P(lu) = K/5Re\v(u;) ■ e*(uj)\ = lu M"(lu) |e(w)| 2 (18) 

which means -P(w) = ^pv 2 \u>\ for a superfluid velocity 
oscillation of amplitude v s . 

V. Numerical Simulation 

Since the predictions of the linearized theory are essen- 
tially classical, we should expect to see the aforemen- 
tioned effects of phonon radiation by solving the non- 
linear Schrodinger equation (NLSE), 



iip 



±VV + (M 2 -1H' 



where we now measure all distances in units of £ and 
times in units of £/c, and ip itself in units of y/ p/m. 
The NLSE was numerically integrated on a two- 
dimensional grid using a operator splitting method |30| . 
To impose the background oscillating superflow the con- 
densate was defined as 



ip{x,t) 



_ iv B (t)-x 



where v s (t) is the chosen time dependent superflow (in 
units of c) uniformly defined over the whole region, and 
cp(x,t) is, initially, a static vortex solution of the NLSE. 
With the initial state representing a condensate with a 
vortex plus a superflow given by v s (t — 0), tp(x,t) was 
evolved according to the NLSE. The resulting equation 
for <p(x, t) is 

up = (ir B • x)<p - i(v + iv s ) 2 <p + (\<f\ 2 - i)<p . (19) 

When v s is constant in time, Galilean invariance means 
that a solution to |l9| is given by 



tp(x, t) = exp(-it/ 2 </2) f(x - v s t) 
0--iV 2 / + (|/| 2 -l)/ . 



In particular, if initially f(x) is a static vortex solution 
of the NLSE, the time evolution of ip(x,t) represents a 
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vortex being rigidly translated with velocity £4 . The devi- 
ations from this "massless" behavior will become evident 
as v s becomes time-dependent. 

In the method used the RHS of |l^ is split into two parts 
which are successively integrated, making the algorithm 
first-order accurate in time [ j30f : 

iCp = — |V <p — iv s ■ V<£ (1st step) 

iCp = v s ■ x if + |u s tp + (\<p\ 2 — l)ip (2nd step) . 

The first step was integrated using the Crank-Nicholson 
method |30| , which is unconditionally stable and second- 
order accurate in time; the second step was integrated ex- 
actly. The time step was 0.01 £/c and the time-dependent 
velocity was along the longest dimension of the grid, a 
channel whose spatial dimensions were 256 x 400 points, 
with a spacing equal to 0.1 or 0.05 £. Along the edges we 
adopted von Neumann boundary conditions for ip(x,t), 
which means that the superfluid velocity computed from 
ip(x,t) was v a (t) at the begining and at the end of the 
channel (this also implies that the initial circular vortex 
flow field had to be slightly distorted). 
There are other ways of imposing the superflow, for in- 
stance using the boundary condition n ■ — iv a ip for 
the full condensate. The influence of the edges however 
takes some time to reach the center where the vortex is 
located, which can be specially inconvenient if one wants 
v s to vary rapidly with time, another problem is that the 
velocity field one obtains is not spatially uniform and 
does not lend itself so easely to a comparison with the 
eletrodynamical theory. 

Fig. 2 shows the effect of a constant acceleration on an 
initially uniform condensate. Here one sees a wake in 
the superfluid density propagating at the speed of sound 
from the beginning of the channel towards the end, and a 
counterwake propagating in the opposite direction. This 
consequence of the imposed accelerated flow takes a finite 
amount of time to reach the channel center, where the 
vortex was placed in the subsequent simulations, and is 
responsible for an observed delay in the vortex response. 
This effect is related to the finite compressibility of the 
superfluid and does no harm to the observation of vortex 
oscillations, it is in fact what causes them. 
A more dangerous effect is the reflection of the wake off 
the end of the channel. To avoid it, one has to restrict the 
observation time to about 40 £/c, for a channel length of 
40 £. This reflection probably explains problems related 
to vortex shedding that were observed in some cases at 
the end of the observation period. 

Results of the Simulations 

The equation of motion one gets for the vortex position in 
the presence of a time dependent background superflow 
is (n = ±1) 



-iujM{u>)V{u>) = ±np z x v s (lo) — V(ui) 



Thus, in the absence of the inertial term, the vortex 
drifts with the superflow. M(u>) is the Fourier trans- 
form of the causal kernel, as above. These dynamics 
imply that for a monochromatic flow one should get a 
response only at the driving frequency or at the reso- 
nance where lo 2 M 2 (uj) — K 2 p 2 , which for a frequency- 
independent mass corresponds to cyclotron oscillations. 
We considered different forms of time dependent flow and 
compared the observed trajectories with the equations of 
motion above. 

To simulate an oscillating flow we took v s (t) to be 



v s (t) 



sin(27rf/T), 



where v° ranged from 0.1 to 0.3 c and T ranged from 2 
to 30 £/c. The trajectories we have obtained display a 
nearly periodic structure with a characteristic frequency 
equal to the driving frequency. We obtained the fre- 
quency dependence of M(ui), as implied by the model's 
equation of motion, by taking the Fourier transform of 
the trajectory and reading its amplitude at the driving 
frequency according to 



M(u>) — ^finp 



uuY{lo) 



(20) 



The results are shown in Fig. 3 together with the func- 
tional form obtained from the linearized effective QED 
theory. 

Since we did not get a perfect steady state response, as 
can be seen in a typical trajectory as in Fig. 4, it was 
not very clear which region of the data array to use in 
taking the Fourier transform. We chose to ignore an ini- 
tial structureless region, corresponding to the delay men- 
tioned above, and to take several Fourier transforms us- 
ing time intervals equal to an integer number of periods, 
all starting at the same point. Hence, with the same tra- 
jectory, we obtained several values of M(lj), which were 
averaged. The error bars in the figures correspond to this 
averaging process. 

We obtained M(w) for 3 amplitudes of the oscillatory ve- 
locity field, namely 0.3, 0.2 and 0.1 c. The last 2 cases 
gave very similar values for the mass, except at lower 
frequencies, where the vortex response was further from 
a steady state, and we could not observe several peri- 
ods of oscillations. Nonetheless, the expected qualitative 
behavior was observed. 

The higher amplitude case, v° = 0.3c, posed more diffi- 
culties at low frequencies. Up to the lowest frequency we 
were able to get, one sees a qualitatively different behav- 
ior, which may be due to the onset of nonlinearities not 
visible in the other cases. 

Note that in the QED theory the real and imaginary parts 
of M(u>) obey a Kramers-Kronig relation. We could not 
check if such a relation existed in the measured M(u>) 
because the points obtained were too scattered to allow 
for a reliable fit. 

We also tried a pulse form for v s (t) , 
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v s (t) = -V2e v max (|;) e 



(21) 



This flow would produce a gaussian displacement along 
the channel for a "massless" vortex. What is observed 
in Fig. 5 is a delayed main peak in the parallel direction 
with an accompanying structure in the perpendicular di- 
rection. The real part M'(oj) is displayed in Fig. 6. Since 
the channel is finite in both width and length, the zero 
frequency limit of M(u>) should cross over to a finite value 
given by the system size. In our system, with length 40 
£, values of M'(u>) greater than In 40 w 3.7 are difficult 
to interpret. 

VI. Impurities and Vortices 

We consider as a simple model of an impurity a point 
object which couples linearly to the boson density, p — 
p(l — B/c). The Lagrangian describing the impurity is 
taken to be 



P 



U a (\x - R a \) B(i 



where R a (t) is the position of the a th impurity. Upon 
integrating out the gauge field A^{x, t), we obtain the 
effective action 

Seff[{*i(t)}, {#»(*')}] = 

fdt I -\kpY, ni e a ? X? If + l m aRl 1 



p K 

Air 



dt dt 



, f d 2 k fduj eW*-*") 



(2ir) 2 J 2ir uj 2 -uj 2 {k) 1 
\p k 2 Y n inj e"* <*' W"^ C» zxfc (t) zxfc-^ (*') 



■— V ni Ifcl C/ a (fc) e -*(^(*)-^(*')) zxk-XAt) 



infrared. The self-interaction term contributes a mass 
shift 



Am a (u) = — J 



rf 2 fc 



fc 4 |/7a(fc)| 5 



m 2 7 (2tt) 2 uj 2 (k) (uj 2 (k)-uj 2 ) 
whose imaginary part (taking uj(k) — cfc), 



(23) 



4m 2 c 6 



\U a (uj/c)\ 2 LU 2 



is super-Ohmic provided U(k — > 0) oc k a with c < | 
§• . 

Consider now the vortex-impurity interaction term. Wc 
will ignore retardation effects, and further assume a point 
interaction between impurities and superfluid, so that 
U a {k) = U a is a constant. Without loss of generality, 
we may consider a single vortex-impurity pair. The con- 
tribution to the effective action is then 



nU pn f i x (X - R) 



2irmc 2 
nUpn 

27T771C 2 



dt 



[X-Rf 
d 



■X 



dtX~0{X - R) 

dx 



where n is the integer charge of the vortex, U is the 
impurity-boson density coupling, and 0(x) = taii~ 1 (y/x) 
is the angle function. Note that if R(t) is time- 
independent, then, for closed paths, S v -i is a topolog- 
ical quantity equal to —(nU p K/mc 2 )W v -i, where W Y -\ 
is the winding number of the vortex about the impurity. 
This makes excellent sense: the vortex effectively counts 
the number of bosons it encircles, and the density re- 
sponse of the superfluid generates a point accumulation 
of AN = —Up/m 2 c 2 bosons to "screen" the impurity. 
Now let the impurity move throughout the superfluid. 
Varying 6> v _i with respect to X a (t) 1 we find the force on 
the vortex to be 



pa _ 



nUpn J 5 af} - 2A a A< 3 



SX C 



2irmc 2 



A 2 



+ ^^fc 2 C4(fc)C4(-fc)e-*^W-^(*'» } (22; 

where ui 2 (k) — c 2 k 2 (l + \^ 2 k 2 ) is the dispersion relation 
derived from the nonlinear Schrodinger equation. (Su- 
perscripts in Eq. |2^ refer only to spatial indices.) We 
now focus on the final two terms, which involve the im- 
purity coordinates. 
The last contribution in Eq. 



22 describes the interac- 



tion between impurities mediated by phonons, including 
a self-interaction term analogous to the self-interaction 
of vortices already discussed. Ignoring retardation ef- 
fects, the purely local superfluid response (see Eq. |^) 
means that the impurity-impurity interaction will be 
short-ranged provided the U a {k) are nonsingular in the 



with A = X — R. This is simply related to the backflow 
current. If we ignore the self-interaction term for the 
vortices, which is appropriate at very low frequencies, 
then the vortex equation of motion, 

npKe^X 13 = F a 

just says that the vortex moves in the dipolar backflow 
field of the moving impurity. 

Finally, consider the force on the impurity due to a mov- 
ing vortex. We find, neglecting the impurity mass renor- 
malization term, 



77 (I 

o-fl 



nUpn \S al3 -2A a A' 3 



2'Kmc 2 



A 2 
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This too has a simple interpretation. The quantity 
Up(R)/m is the local potential in which the impurity 
moves, and hence the force on the impurity is F = 
— (U/m)\/p. Now recall Duan's result (Eq. |]) for the 
density response to a moving vortex. Taking the gradi- 
ent gives us the appropriate force. 

Polaron Model of a Quantum Vortex 

Niu, Ao, and Thouless [NAT] have investigated a model 
of a quantum vortex coupled to superfluid density fluctu- 
ations ||. They describe the vortex as a nonrelativistic 
particle of mass M c in a background uniform magnetic 
field (corresponding to the average superfluid density) 
and assume a scalar coupling to the phonons, i.e. 

Hnat = 2M ^ + nK P d / c "> 2 + E huJ k( d l d k + 

k 



-Y,W(k)e^(d z + dl,) , 



where n is the integer vorticity, Vxa = — cz accounts for 
the geometric phase due to the background bosons, and 
ujj. is the phonon frequency at wavevector k. NAT go 
on to investigate a simple polaronic wavefunction which 
accounts for the phonon cloud around a vortex and con- 
clude that an infinite vortex mass would shrink the quan- 
tum uncertainty in the vortex position to zero, a situa- 
tion in conflict with explicit calculations using Feynman's 
trial vortex wavefunction. At this level, the mass M e 
is phenomenological and does not include the effects of 
phonons. It may be perhaps more appropriate to con- 
sider M e as the mass of an external particle trap ped in 
the vortex core, as considered by Demircan et al. |3l[| . 
In keeping with the general philosophy that the vortex is 
a topological object which couples to the boson density 
as in Eq. H, we propose a variant of the NAT model: 



H= 2M~ (p + nK P( S +A)/c) +J2 



where u) r = Tick and 

k 



fc v k k 2 y ' 



/ tic 
2p\%\ 



x k 



ik- 



d lk) 



is the quantized radiation field corresponding to the su- 
perfluid density and current fluctuations p6[ (f2 is the 
area of the system). We rewrite TC as 



n = 



n 2 

2M r 



nnp 



An 



1 ? —? 

n k p 
2M c c 2 



vortex of strength n. We now work to lowest order in A, 
and follow NAT by assuming a trial state 



with 



,-(f-S.f/il 3 -iz-rxR/2l 2 a 



treating I as a variational parameter. Taking the expecta- 
tion value of Tt is the state X^(r), we obtain the effective 
phonon Hamiltonian 

H & = E K-KA- + 1) + E w &) ( d k + d -^ 



k 



phujr 1 



The phonon ground state is a coherent state, 
W(k) 



k k 



and the total energy is 



E : 



2M 



k 



| + (l-nVM e )^ 



where p — irpS, 2 — k 2 p/Attc 2 as before. Setting dE/dl = 
gives 



Note that no solution exists for M e < n 2 p, which we 
interpret in the following manner. A cyclotron mode will 
show up as a pole in the denominator of Eq. [l7|, which 
means ujM c (lo) = ±Kp~. Now no such pole exists in the 
absence of the external mass M e , but if M c is added 
to our M(oj), then a damped cyclotron resonance does 
exist at the cyclotron frequency uj ~ np/M c , provided 
that M e > n 2 p. 

NAT compute a renormalized magnetic length I accord- 
ing to the relation 



where 77 = p — ^nnpz x r is the cyclotron momentum 
operator for the vortex. Note that [n x ,n y ] = ih 2 /l 2 , 
where l a = y/m/2wnp is the "magnetic length" for a and 



>(.R) \V(R + f))) =exp 

which is obtained from the overlaps 

iz ■ R x fj 



-m-+o(\v\ 4 ) 

2l 2 



exp 



211 



8l 2 
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exp — > J — ( 1 — cos k ■ f) 



We find 



irn ( m \ 2 £ 3 Z 



1 



i i 

P ~ 2P + 
VII. Conclusion 

In this paper we have explored the theory of dynam- 
ical vortices in superfluid films, deriving a frequency- 
dependent vortex mass which enters into the vortex equa- 
tions of motion, as well as describing dissipation by ra- 
diation of phonons. Numerical simulations corroborat- 
ing the predicted behavior of M(u) were presented as 
well. These calculations may be extended to (3 + 1)- 
dimensional superfluids as well These results will 

be presented in a future publication. 
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Appendix: Inertial Term for Many Vortices 

The velocity-dependent part of the effective action in- 
duced by integrating out the phonon field is 



AS = -\pn 2 dt dt' 



d 2 k f did 



j(t-i') 



niUj e 



(2tt)V 2vr lu 2 - LU 2 (k) 



-i&iXiW-Xjit')) ix j,. x ^ t ) i x k.x 3 (t') . 



We define A tj {t, t') = X t (t) - X^t'). Using 
dk 



2tt 



e- lk ' A k a k p = \j (kA)5 al3 



+ U 2 (kA)(S af3 -2A a A 13 ) 



where J n {z) is the Bessel function of order n. Taking 
u>(k) = ck, we have 



AS = - 



dt dt' =- e""- 1 '' Yl n ^ (*') 



8ttc 2 J J J 2n 
K 2 {-iuj Ai 3 1 c) + 



2c 2 



uj 2 A 2 - 



(5» u ~ 2A»A V ) 



Koi-iuAij/cjS^ 



where K n (z) is a modified Bessel function, and where 
lu — > lu + i0 + is understood. Since the above integrand 
is already quadratic in velocities, which we assume are 
small compared with c, we may approximate Aij as a 



constant. Expansions of K n {z) for small z 25 yield 

K Q (-iz) = -C - hx(-iz/2) + ... 
2 

K 2 (-iz) + ^ = \z 2 \n{-iz/2) + ... 

and at low frequencies the first of these terms dominates, 
so provided lu <C c/A lms we recover the action of Eq. fill 
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FIGURE CAPTIONS 



FIG.l Complex frequency-dependent vortex mass 
M(ui) from eq. |l4|. Also shown is the complex dimen- 
sionless inertial parameter r(w). 

FIG. 2 Amplitude of 4>{x,t) for a uniformly accelerated 
condensate, plotted as a function of length along the di- 
rection of flow, at ten equally spaced time intervals sep- 
arated by At = 5£/c. 

FIG. 3 Complex mass M(ui) inferred from eq. |o[ Solid 
curves are M'(oj) and M"(lo) from the linearized electro- 
dynamic theory, (a) M'(w) for driving amplitude 0.1c 
(b) M"(u) for driving amplitude 0.1c. (c) M'(u>) for 
driving amplitude 0.2 c. (d) M"(lv) for driving ampli- 
tude 0.2 c. (e) M'(uj) for driving amplitude 0.3 c. (f) 
M"{lo) for driving amplitude 0.3 c. 

FIG. 4 Typical vortex trajectory X(t) with harmonic 
forcing of period T = 5£/c. 

FIG. 5 Vortex trajectory X(t) for a pulsed superfiow 



given by eq. 21. In this case T = 2£/c and v max = 0.4c. 
The smooth curve is the pulse shape. 

FIG. 6 Real part of the mass M'(lo) obtained from anal- 
ysis of response to the pulse flow. Solid curve is the pre- 
diction of the linearized electrodynamic theory. 
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Vortex Mass: Results of Linearized Effective Action 
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